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ABSTRACT 

We have performed smoothed particle hydrodynamics (SPH) simulations to 
study the response of the central kiloparsec region of a gaseous disk to the imposi- 
tion of nonaxisymmetric bar potentials. The model galaxies are composed of the 
three axisymmetric components (halo, disk, and bulge) and a non-axisymmetric 
bar. These components are assumed to be invariant in time in the frame coro- 
tating with the bar. The potential of spherical 7-models of Dehnen is adopted 
for the bulge component whose density varies as r~ 7 near the center and r~ 4 
at larger radii and hence, possesses a central density core for 7 = and cusps 
for 7 > 0. Since the central mass concentration of the model galaxies increases 
with the cusp parameter 7, we have examined here the effect of the central mass 
concentration by varying the cusp parameter 7 on the mechanism responsible 
for the formation of the symmetric two-armed nuclear spirals in barred galaxies. 
Our simulations show that the symmetric two-armed nuclear spirals are formed 
by hydrodynamic spiral shocks driven by the gravitational torque of the bar for 
the models with 7 = and 0.5. On the other hand, the symmetric two-armed 
nuclear spirals in the models with 7 = 1 and 1.5 are explained by gas density 
waves. Thus, we conclude that the mechanism responsible for the formation of 
the symmetric two-armed nuclear spirals in barred galaxies changes from the hy- 
drodynamic shocks to the gas density waves when the central mass concentration 
increases from 7 = to 1.5. 
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INTRODUCTION 



High resolution visible and near- infrared observations by the Hubble Space Telescope 
(HST) and ground-bas ed telescopes with adapti v e optics show nuclear spirals in the cen- 



ters o 



2002: 



Martini et al 



2003a 



as well as in those of normal galaxies (IPhillips et a l. 1996; 



sm 

active galax i es (jRegan fc Mulchaeyl Il999l ; iMartini k, Poggd 11999c IPogge fc Martini 



Elmegreen et alill998l : lLaine et alill999l . l200ll ; ICarollo et al1l2002l ; IMartini et al.l l2003al lbh. 



Nuclear spirals have a variety of morphologi es such as grand-design spirals, one-armed spirals, 



and flocculent spirals (IMartini et al.ll2003al ). The grand-design nuclear spirals are symmetric 
two-armed spirals, while the flocculent, or "chaotic," nuclear spirals have multiple arms that 
are likely to be segmented into smaller arms. 

The morphology of symmetric two-armed nuclear spirals depends on the gas sound speed 

and the shape of the gravita tional potential in the nuclear regions of gal axies (lEnglmaier fc Shlosman 
2000; lAnn fc Thakur 2005 ). A number of hydrodynamic simulatio ns ( Englmaier fc Gerhard 



199/ 



Patsis &: Athanassoulal I200CH : lEnglmaier fc Shlosman! |2000| ; iMaciejewski et al. 



2002 



Ann fc Thakur 20051 ) have already been performed to understand the effect of the gas sound 
speed on the gas flows inside the in ner Lindblad resonance (ILR ). Th e smoothed particle hy- 
drody namics (SPH) simulations of lEnglmaier fc Gerhard! (119971 ) and lPatsis fc Athanassoula 
(120001 ) showed that the off-axis shocks form at low gas sound speed, while the on-axis shocks 
develop at high gas sound speed. However, they were unable to analyze the nuclear mor- 
phology in detail due to the lack of resolution of their simulations. On the other hand, the 
symmetric two-armed nuclear spirals generated by the grav itational torque of the bar wer e 
found in the high - resolu tion hydrodynamic simula tions of Englmaier fc Shlosman! (2000), 
Macieiewski et al.l ( 120021 ) , and I Ann fc Thakur! (120051 ) . It was shown by lEnglmaier &: Shlosman 
(120001 ) with grid-based hydrodynamic simulations that the symmetric two-armed nuclear spi- 
ral can be explained in terms the gas density waves, since the spiral arms open out as the 
gas sound speed increases. A similar effect of the gas sound speed on the o penness of sym- 
metric two-armed nuclear spirals was also found in the SPH simulations of lAnn fc Thakur 
(120051 ) . However, their simulations suggested that the symmetric two- ar med nuclear spirals 



are for med by hydrodynamic spiral shocks driven by the bar. Furthermore, IMaciejewski et al. 



(120021 ) found that the gas sound speed plays a crucial role in the mechanism responsible for 
the formation of the symmetric two-armed nuclear spirals in barred galaxies. Their grid- 
based hydrodynamic simulations showed that the symmetric two-armed nuclear spirals are 
formed by hydrodynamic shocks driven by the bar when the gas sound speed is high (c s = 20 
km s _1 ), while for gas at a low sound speed (c s = 5 km s" 1 ), they are mai ntained by g a s den- 
sity waves. It was noticed in the grid-based hydrodynamic simulations of lMaciejewskil (120041 ) 
that like the gas sound speed, the bar strength can also affect the mechanism responsible for 
the formation of the symmetric two-armed nuclear spirals in barred galaxies. 
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Since the central mass concentration determines the shape of the gravitational potential 
in the nuclear regions of galaxies, it can also affect the morphology of symmetric two-armed 
nuclea r spi rals. This c an be understood from the simulations of lEnglmaier fc Shlosman 
( I2000l ) and lAnn &: Led (120041 ). where they have mentioned that the openness of symmetric 
two-armed nuclear spiral arms d epends on the central mass concentration. Our previous 
simulations (lAnn fc Thakurll2005l ) showed that the presence of a central supermassive black 
hole (SMBH), which can greatly affect the central mass concentration, plays a significant 
role in the sense that the symmetric two-armed nuclear spirals form when the mass of an 
SM BH is massive enoug h to remove the inner inner Lindblad resonance (IILR). In addition to 
this, iMaciejewskil (120041 ) noticed that a central massive black hole (MBH) of sufficient mass 
can allow the symmetric two-armed nuclear spirals to extend all the way to its immediate 
vicinity. 

However, it is important to note here that the effect of the central mass concentration 
on the mechanism responsible for the formation of the symmetric two-armed nuclear spiral 
was not known from pr evious studies, whereas t hat of the gas sound speed was understood 
from the simulations of iMaciejewski et al.l (120021 ) . Furthermore, it is worth mentioning that 
if we consider the model galaxy whose bulge central density varies as r -7 , its mass would 
be more centrally concentrated as the cusp parameter 7 increases. Thus, the purpose of 
the present paper is to examine the effect of the central mass concentration by varying the 
cusp parameter 7 on the mechanism responsible for the formation of the symmetric two- 
armed nuclear spirals in barred galaxies. For this purpose, we employ the SPH technique to 
solve the hydrodynamic equations for the gas flow in the model galaxies whose bulge central 
densities vary as r~ 7 . We have also included the self-gravity of the gas in our calculations, 
since this can play a significant role in high-density regions such as shocks. 

The remainder of this paper is organized as follows. In § 2 we describe the model galaxy, 
model parameters, the gaseous disk model, and the numerical method. In § 3 we present the 
results of our simulations, and the discussion of our results is given in § 4. The conclusions 
are presented in § 5. 



2. MODELS AND NUMERICAL METHODS 
2.1. The Model Galaxy 



The model galaxy is assumed to be composed of four components: halo, disk, bulge, and 
bar. The properties (density, size and structure) of all these potential generating components 
are assumed to be invariant in time in the frame corotating with the bar. 
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The halo component, which gives rise to the flat rotation curve at outer radii, is assumed 
to have a logarithmic potential with finite core radius, 

1 



halo 



-v \n(R h + r ) + const, 



1 



where Rh is the halo core radius and Vq is the constant rotation velocity at large r. 



The disk component is assumed to have the potential of Freeman's! ( 119701 ) exponential 



disk, 



$(r) 



disk 



- h 



2R n 



2R« 



2R n 



2R„ 



(2) 



where R^ is the disk scale length, E is the central surface density, and Iq, K , Ii, and K\ 
are the modified Bessel functions. The total disk mass is simply M^sk — 2vrS -R^. 



For the bulge component, we adopt the potential of the spherical 7-models of iDehnen 



(1993) 



GM, 



bulge 



bulge 



(2 - 7)^6 



r + R t 



1 



, for 7 ^2, 



(3) 



where M^uige is the total mass of the bulge, Rb is the scale-length, and 7 is the cusp parameter 
that can have a value in between < 7 < 3. Its corresponding density varies as r -7 near the 
center and r~ 4 at large radii, and hence possesses a central density core for 7 = and cusps 
for 7 > 0. Th is bulge with 7 = 1.5 most closely resembles the de Vaucouleurs' i? 1//4 -profile 
( )Dehnenlll993l ) that is known to better represent the bulges of real galaxies. Hereafter, the 
central mass concentration will be represented by 7, since the central mass concentration of 
the model galaxies increases with the cusp parameter 7. 

The bar is a triaxial component in three dimensions. However, our simulation is re- 
stricted to the two-di mensional disk. Thus, w e assume the following form of the gravitational 
potential proposed by lLong &: Muralil (119921 ) for the bar component: 



$(r) 



GM, 



bar 



bar 



2a 



log 



x — a + T_ 
x + a + T + 



(4) 



where T± = \J[(a ± x) 2 + y 2 + 6 2 ], and a and b are the needle length and softening length, 
respectively. The parameter a/b defines the elongation of the bar, which is propo rtional 
to the axis ratio a/b of the flattened hom ogeneous ellipsoid bar of Freeman (jl966 ) when 
a/b > 2 (see Fig. 3 of lLong &: Muralil Il992l ). where a and b are the major and minor axes of 
a flattened homogeneous ellipsoid bar. 
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TABLE 1 



MODEL PARAMETERS 





c s 




Model 


(km s _1 ) 


7 


Ml 


10. 


0.0 


M2 


10. 


0.5 


M3 


10. 


1.0 


M4 


10. 


1.5 


M5 


15. 


0.0 


M6 


15. 


0.5 


M7 


15. 


1.0 


M8 


20. 


0.0 


M9 


20. 


0.5 


M10 


20. 


1.0 
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We have considered ten models M1-M10, which are listed in Table 1. All these models 
are defined according to different choices of the gas sound speed c s and the central mass 
concentration 7. In this paper, we are mainly interested to examine the effect of the central 
mass concentration 7 on the mechanism responsible for the formation of the symmetric two- 
armed nuclear spirals in barred galaxies when the gas sound speed is chosen to be c s = 10 km 
s _1 (see § 3.2 and § 3.3 ). Thus, apart from the models with 7 = 0, 0.5, and 1, an additional 
model with 7 = 1.5 is considered for the gas sound speed of c s = 10 km s -1 (i.e., model 
M4). However, the models with 7 = 0, 0.5, and 1 are chosen for each of the remaining gas 
sound speeds of c s = 15 and 20 km s . Since the total visible mass of the model galaxies is 
assumed to be Mq ~ 4x 10 10 , the masses of the gas, disk, bulge, and bar components in units 
of M sc = 2x 10 n M Q are M gas = 0.002, M disk = 0.104, M bulge = 0.054, and M bar = 0.04, 
respectively. Thus, the mass distributions of the model galaxies are chosen to be similar 
to those of early- type barred galaxies (~SBa), since the disk-to-bulge mass ratio and the 
fractional mass of the bar are Mdisk/^buige = 2 and M bar /MG = 0.2, respectively. The scale 
lengths of the halo, disk, and bulge are assumed to be Rh = 15.0 kpc, Rj, = 3.0 kpc, and 
Rb = 0.25 kpc, respectively. Furthermore, the bar needle length a, the bar softening length 
b, the bar pattern speed Q p , and the rotation velocity of the halo Vq are fixed at 3.0 kpc, 1.0 
kpc, 42.8 km s _1 kpc -1 , and 186 km s _1 , respectively. 

Figure 1 shows the rotation curves and angular frequency curves of the models mentioned 
in Table 1 that have the central mass concentrations of 7 = and 1. Since the models with 
constant values of 7 have the same mass distribution, we present the rotational curves and 
angular frequency curves of models Ml and M3 that assume 7 = and 1, respectively. The 
top left panel of Figure 1 represents the rotational velocities of the halo, disk, bulge, and bar 
components of model Ml with 7 = as a function of radius, while the top right panel shows 
the rotational curves of model M3 with 7=1. Since the bar is nonaxisymmetric component, 
i ts contribution to the rotational velocity is included here after averaging it axisymmetrically 



( iHasan et al.lll993l ). As can be seen in the top panels of Figure 1, the rotational velocities in 
the vicinity of the galactic center are affected by the increase in central mass concentration 
from 7 = to 1. However, the bottom panels of Figure 1 represent the angular frequencies 
corresponding to the rotational velocity generated by the axisymmetric components as a 
function of radius. The adaption of Q p = 42.8 km s _1 kpc -1 places the IILR at the negligible 
distance from the galactic center, and the outer inner Lindblad resonance (OILR) near 1.3 
kpc in model Ml with 7 = 0, whereas it places the single ILR close to 1.5 kpc in model M3 
with 7 = 1. Furthermore, the corotation radius (Rcr) is found to be around 1.2a in models 
Ml and M3 for the above choice of Q p . Since the angular frequency curve reaches high values 
(larger than Q p ) at the nucleus for model Ml that considers the central mass concentration 
of 7 = 0, it suggests that for the assumed values of Mb u i ge = 0.054 and Rb = 0.25 kpc, the 
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mass of model Ml with 7 = is sufficiently centrally concentrated (IMaciejewskil 120031 ) to 
weaken the IILR by placing its location at the negligible distance from the galactic center. 
However, the IILR is completely removed and the angular frequency curve goes to infinity at 
the galactic center for model M3 with 7 = 1, since the central mass concentration increases 
from 7 = to 1. Although it is not presented here, we would like to mention that except for 
the slight difference in the very vicinity of galactic center, the rotation curves of the models 
with the remaining values of 7 = 0.5 and 1.5 appear to be similar to those of model M3 with 
7 = 1. In addition to this, the angular frequency curves of the models with 7 = 0.5 and 1.5 
do not greatly differ from those of model M3 with 7 = 1, except for the minor shift in the 
location of the single ILR. 



2.3. The Gaseous Disk Model 



In the direction of modeling the response of the gaseous disk to the imposition of 



nonaxisymmetric bar pote ntials, many significant resu 



ts have been produced so far when 



the gas is isothermal (eg., 


Englmaier & Gerhard 


1997: Fukuda et al. 


1998. 


2000; 


Fux 


1999. 


2001; 


Lee et al. 19991; Patsis & Athanassoula 


2000|; 


Enelmaier & Shlosman 


2000; 


Ann 


2001; 


Mack 


newski et al. 2002; iMacieiewski 


2003. 


20oi 


Ann & Lee 2004: 


Ann & Thakur 


2005; 


Perez 


2008: 


Namekata et al. 2008). 


Endmaier & Gerhard! (1997h. 1 


3nelmaier & Shlosman 


(2000 


), and 


Maciejewski (! 


2003) pointed out that this approach of considering isothermal gas 



is consistent with the ISM model of ICoweil (119801 ) . who argued that the cloud fluid can be 



treated like an isothermal gas if the clouds have an equilibrium mass spectrum. The equilib- 
rium is assumed to be maintained by a steady supply of small clouds by supernovae. Thus, 
we assume the infinitesimally thin and isothermal uniform gaseous disks with the effective 
sound speeds of c s = 10, 15, and 20 km s -1 . The mass of the gaseous disk was chosen to be 
1% of the total visible mass of the model galaxy, since we simulate the response of the central 
kiloparsec region of the gaseous disk in early-type barred galaxies. In order to represent the 
initial gaseous disk with constant density, 2 x 10 4 SPH particles are distributed uniformly 
over a grid inside a circular disk of 5 kpc radius. The initial circular velocities are given 
to the SPH particles so that the centrifugal accelerations of the particles are balanced with 
their gravitational accelerations. In the context of this equilibrium gaseous disk model, the 
pressure gradient forces are zero. The initial resolution length, defined by equation (12), 
was found to be ~ 147 pc. However, the resolution of the SPH simulations becomes far 
better than the initial values in the central region as the density in the nuclear region in- 
creases within a few bar rotation periods (see § 2.4). We introduce the bar potential slowly 
within one-half bar rotation period to avoid a spurious response of the gaseous disk, where 
the bar rotation period r^ ar is ~ 1.44 x 10 8 yr. The previous simulations of the isother- 
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mal gas flow in barred galaxies with broad spectrum of values for effective sou nd speed 



( lEnglmaier &: Gerhard! 1 19971 ; iPatsis fc Athanassoula 



200(1 iMacieiewski et al.l 12002( 1 and the 



compilation of HI velocity dispersion by iKamphuid (119931 ) suggest that our adopted values 
of c s = 10, 15, and 20 km s _1 correspond to the vertical cloud velocity dispersion in the 
Galactic bulge. 



2.4. Numerical Method 



We have used the PMSPH method to solve the hydrodynamical equations for the gas 
flow in barred galaxies. Our simulations are confined in two dimensions, since restricting 
the calculation to two dimensions increases gr eatly the resolution, which is important in the 
present problem of re s olving the shocks (cf. lEnglmaier fc Gerhard! 1 199 71 ). In this regard, 
Englmaier fc Gerhard! (119971 ) mentioned that the two-dimensional description appears to be 
a reasonable approximation, since they have done some three-dimensional simulations, where 
the gas flows were inde ed similar to the two-dimensional ones. We h ave used the PMSPH 
code of lFuxl (119991 . l200ll ). which has recently been used by iPerea (120081 ) to model the gas flow 
in the barred galaxy NGC 4123. R. Fux of Geneva Observatory Group has kindly provided 
the PMSPH code to us. Here we briefly describe the necessary sp ecifications relevant to 
the particular aspects of our PMSPH code and re f er the read e rs to iFuxl (J1999, 200l|) for a 



detailed description of our code and to lBenzl (Il990l ) . iMonaghanl (119921 ). and lSteinmetzl (119961 ) 



for the general properties of solving Euler's equation of motion using the popular Lagrangian 
method. 



2.4-1- Gas Hydrodynamics with SPH 



Since SPH is fully Lagrangian, it does not have any grid to impose any artificial re- 
strictions on the global geometry of the systems under study or any mesh-related limitations 
on the dynamic range in spatial resolution. Thus, a high nume rical resolution is n aturally 
achieved in high-density regions such as shocks by using SPH (IFukuda et al.lll998l ). which 
is necessary to solve our present problem of understanding the response of gaseous disks in 
the nuclear regions of barred galaxies. 

Using the SPH method, we solve the Euler equation of motion that can be written as 

v(p + n) 



dv dv , „. 

* S 9F + (V ' V)V 



(5) 



where v(r) is the velocity field and p(r), P(r) and n(r) are the density, pressure, and artificial 
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viscosity of the gas, respectively. — V$ is the gravitational force, where $ includes potentials 
of the g well as the stellar and dark components. 

In SPH, the density is evaluated directly at each particle position in space by summing 
the contributions from the density profiles of neighboring particles and is represented by 

Kg 

Pi = p(Ti) = ^mjWfa - rj, h), (6) 

i=i 

where N g is the number of gas particles, rrij is the mass of individual particles, W(r, h) is 
the kernel function, and h is the smoothing length that defines the local spatial resolution 
and is proportional to the local mean interparticle spacing. In our code, the adopted SPH 
kernel is a spherically symmetric spline that vanishes outside 2h because of its finite spatial 
extension. 

Our code assumes an isothermal equation of state, i.e., 

Pi = <? s Pi, (7) 

where c s is the effective sound speed, which can be interpreted globally as the velocity 
dispersion of the interstellar clouds. The internal energy of the gas is at any time and 
everywhere constant, and in particular, there is no energy eq uation involv ed in our code. 



With the above equation (7), the pressure gradient term of iBenzl (jl990l ) in the Euler's 
equation simply takes the form 



(?) 



VV- - + — ViW{ri - rj, h). 

VP* Pjj 



The introduction of an artificial viscosity is required for an accurate tre atmen t of th e 
flow near shocks. In our code, the artificial viscosity is calculated exactly as in IBenzl (Il990l ). 
It consists of the bulk and the von Neumann-Richtmyer vis cosities, w i th the standard pa- 
rameters set to a = 1.0 and (3 = 2.5, and takes into account iBalsara'd (119951 ) correction to 
avoid energy dissipation in pure shearing flows. 

In order to achieve high resolution, it is necessary to treat the smoothing length h 
properly, since it determines the local spatial resolution. To increase the resolution in high- 
density regions like shocks, our code uses a spatially variable smoothing length h to assign 
an individual smoothing length hi to each particle in such a manner that the number of 
neighbor particles iVj to each particle always remains as close as possible to a fixed number 
N D . This allows us to save significant computing time while keeping spatial resolutions high 
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enough to resolve the shocks (lFuxlll999l . 120011 ). Two particles i and j are defined as mutual 
neighbors if i 7^ j and |rj — Tj \ < 2hij, where hij = (hi + hj) /2 is the symmetrized smoothing 
length, which is supplied in equations (6) and (8), as well as in the formula for the artificial 
viscosity to ensure momentum conservation. 



At each time step, the smoothing lengths are updated in our code according to the 
general three-dimensional scaling law: 



hi 

hn 



Nj + l pc 
N + l Pl 



1/3 



(9) 



where h Q and p a are constants, iV$ is th e number of n eighbors of particle i, and +1 is added to 
take into account particle i. Following iBenzl (119901 ) . we take the time derivative of equation 
(9) and substitute the continuity equation to yield 



• dhi 1 
hi = — — = -hi 
dt 3 



dNi 



Ni + 1 dt 



+ [V 



(10) 



In equation (10). iBena (119901 ) did not include the term dNi/dt to ensure a constant number 
of neighbors. However, this does not prevent a slow numerical departure of the iVj's from N a 
with time. To overcome this problem, our code simply damps such departures by setting 



dNi N - N n 



dt 



r]At 



in equation (10) and integrates the resulting equation along with the equations of motion. 
The parameter rj controls the damping rate per time step At and should be significantly 
greater than 1 to avoid abrupt discontinuities in nongravitational forces. In all our simula- 
tions, we have assumed iV = 35 and rj = 5. 



The initial smoothing lengths of the gas particles, confined near the plane z 
defined as 



hi(t = t Q ) 



N mj 
4nZ(R 9 ) 



0, are 



(12) 



where m; and £(-R 9 ) are the mass of each gas particle and the surface density of the initial 
gaseous disk of radius R g = 5 kpc, respectively. 

It is customary to impose quite a large lower limit on the minimum smoothing length 
hmin to save computing time. However, because of the fast neighbor searching algorithm 
adopted in our code, we have used h min = 0.1 pc to resolve structures down to the parsec 
scale. Although the initial smoothing length (i.e., the initial resolution length) of ~ 147 pc 
is employed in our simulations (see § 2.3), the smallest smoothing length actually achieved 
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after a few bar rotation periods is ~ 2 pc, which i s sma ller than the h m i n = 5 pc used in the 
high-resolution SPH simulations of iFukuda et al.l (120001 ) . We find that the smoothing length 
of ~ 2 pc is small enough to resolve the nuclear spirals, as shown in § 3. Furthermore, it is 
worth mentioning he re that our code uses a synchronized version of the standard leap-frog 
time integrator (e.g.. lHut et al.lll995l ). which is known to well conserve the total energy and 
the total angular momentum and also employs an adaptive tim e step to temporally resolve 
the high-density shocks in the gaseous disk (see IFuxn1999i . l200fl ). 



2.4.2. Gravitation with PM 



Our gaseous disk is not a self-gravitating one, since its mass is assumed to be only a tiny 
fraction of the model galaxy (~ 1%). This suggests that Toomre' s Q value is much greater 
than unity (i.e., Q >> 1), which represents a stable ga seous disk (cf . lElmegreen fc Thomasson 
19931 : IFukuda et all 1 19981 . I2OO0I : knn fc Thakuri I2OO0} ) . However, we have included the self- 
gravity of gas in our calculations because it can play a significant role in high-density regions 
such as shocks that develop in the later evolution of gaseous disk. To compute to the self- 
gravity of gas, our co de uses the particle-mesh ( PM) technique with polar-cylindrical grid 
geometry described in lPfenniger Sz Friedli Jl993h . The advantage of considering the polar- 
cylindrical grid is that the radial and azimuthal resolutions increase toward the center, where 
a variety of features are likely to develop because of the gas inflow driven by the bar. Using 
the convolution theorem, the potential is computed on the polar-cylindrical grid by the fast 
Fourier transform (FFT) technique in the <fi and z dimensions, where the cells are equally 
spaced. The radial spacing of the cells is logarithmic with a linear core to avoid an accumu- 
lation point at the center. Since the Keplerian kernel does not truly satisfy the conditions of 
the convolution theorem, a variable homog eneous ellipsoidal kernel is adopted in our c ode for 
the softening of the short range forces (cf. iPfenniger &: Friedlil Il993l ; lFuxill999l . l200ll ). This 
kernel tends toward the Keplerian ker nel at large distance, and can also be useful in situa- 
tions that need anisotropic resolution (IPfenniger &: Friedlilll993l ). Since the polar-cylindrical 
grid is non- homogeneous, a variable softening length is necessary. Otherwise a softening 
length intermediate between the smallest and largest cell sizes both decreases resolution in 
some parts of the grid and enhances grid noise in other parts of the grid. Thus, the softening 
lengths represented by the semi-axes of our adopted variable homogeneous ellipsoidal kernel 
have been set as 1.1 times the respective cell dimensions. This ensures the conservation 
of the tota l energy, as well as of the total angular momentum while keeping high central 
resolution (IPfenniger &: Friedh1ll993l ; lFuxill999l . l200ll ). For a mor e detailed description of the 
PM technique used i n our code, we simply refer the readers to IPfenniger &: Friedlil (119931 ) 
and lFuxl Jl999l . l200lh . 
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3. RESULTS 

3.1. Nuclear Feature of the Gaseous Disk 

The nuclear features of the gaseous disks in models M1-M10 at the evolution time of 
10Tb ar are shown in Figure 2. Here our results are shown on the frame corotating with the 
bar, which always lies horizontally. As can be seen in Figure 2, models M1-M10 clearly show 
the development of symmetric two-armed nuclear spirals at the evolution time of 10r;, ar . For 
models M1-M4 in the first column that consider the gas sound speed of c s = 10 km s _1 , 
the innermost parts of the symmetric two-armed nuclear spiral of model Ml with 7 = 
extend up to winding angle of ~ jn, whereas they can wind up to ~ 2n for model M4 with 
7 = 1.5. It is also apparent that the nuclear disks are completely filled with the gas particles 
for models M8-M10 in the last column that assume the gas sound speed of c s = 20 km s _1 . 
In this case, the innermost parts of the symmetric two-armed nuclear spiral of model M8 
with 7 = extend all the way to the center with winding angle of ~ 37r, whereas they appear 
to be tightly wound in the sense that their arms cannot be resolved after ~ |7r for model 
M10 with 7 = 1. These characteristics of the nuclear features show that the symmetric 
two-armed nuclear spiral arms open out as the central mass concentration 7 decreases while 
keeping the gas sound speed c s constant. By comparing the winding angles corresponding to 
the symmetric two-armed nuclear spirals of the models with constant values of 7 in Figure 
2 such as those mentioned above for models Ml and M8 that assume 7 = but different 
gas sound speeds of c s = 10 and 20 km s , respectively, it is clear that the nuclear spiral 
arms also open out as the gas sound speed c s increases while keeping the central mass 
concentration 7 constant. In addition to the above nuclear features of the gaseous disks in 
models M1-M10, Figure 2 also shows that the diameter of the nuclear disk depends on the 
central mass concentration 7 and the gas sound speed c s . The models with high central mass 
concentrations allow a larger nuclear disk than those with low central mass concentrations. 
Furthermore, a larger nuclear disk is also found in the models with low gas sound speeds. 
As a result, a largest nuclear disk among our selected models is noticed in model M4, which 
assumes 7 = 1.5 and c s = 10 km s _1 . 



3.2. Nuclear Distribution of the Shocked Gaseous Particles 



As shown in the previous section, all the models in the present study show the formation 
of the symmetric two-armed nucle ar spirals at the evolution time of 10rfo rtr . The mechanisms 
such as the hy drodynamic shocks (IMacieiewski et al.ll2002l ; lAnn fc Thakurll2005l ) and the gas 
density waves (jEnglmaier fc Shlosmanl|2000l ) were proposed for the formation of such nuclear 
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spirals in barred galaxies. Since a significant role is played by the gas sound speed c s to decide 
whether the symmetri c two-armed nuclear spir als are formed by hydrodynamic shocks or by 
gas density waves (cf. iMaciejewski et al.ll2002l ). we are interested here to examine any such 
role of the central mass concentration 7. For this study, we mainly consider models M1-M4, 
which have c s = 10 km s _1 but different central mass concentrations of 7 = 0, 0.5, 1, and 
1.5, respectively. In order to see whether the symmetric two-armed nuclear spirals shown in 
Figure 2 for models M1-M4 are made of shocked gas particles, we exa mine the distributions 
of ga s particles that are sho cked by supersonic flow, i.e., —hV-v > c s (jEnglmaier fc Gerhard 



19971 ; lAnn &: Thakurl 120051 ) . in the central kiloparsec of models M1-M4 at the evolution time 



of 10Tb ar in Figure 3. As can be seen in Figure 3, the distributions of shocked gas particles 
in the central kiloparsec of models Ml and M2 reveal well-defined spiral patterns whose 
innermost parts extend up to winding angle of ~ 2tt. These nuclear spiral shocks look very 
similar to the symmetric two-armed nuclear spiral patterns developed in Figure 2 for the 
distributions of all the gas particles in the central kiloparsec of models Ml and M2. On the 
other hand, it is also apparent from Figure 3 that the distributions of shocked gas particles 
in the central kiloparsec of models M3 and M4 cannot move inward in a curling manner after 
winding angle of 71. However, the distributions of all the gas particles in the central kiloparsec 
of models M3 and M4 show the development of symmetric two-armed nuclear spirals whose 
innermost parts extend up to winding angle of ~ 2tt, as shown in Figure 2. Thus, our 
simulations suggest that the symmetric two-armed nuclear spirals of models Ml and M2 
that respectively assume 7 = and 0.5 are formed by hydrodynamic spiral shocks driven by 
the gravitational torque of the bar, whereas they cannot be supported by the shocked gas 
particles for models M3 and M4 that respectively assume 7 = 1 and 1.5. It is worthwhile 
to mention here that the models with each of the remaining values of the gas sound speed 
(i.e., c s = 15 and 20 km s _1 ) show a similar effect of the central mass concentration 7 on the 
formation of the nuclear spirals in barred galaxies to that shown in Figure 3 for the models 
with c s = 10 km s . 



3.3. Effect of Central Mass Concentration on Gas Inflow 



Symmetric two-armed nuclear spirals are suggested to be formed by the gas inflow across 
the ILR du e to the loss of angular momentum at the shoc ks driven by the gravitational torque 
of the bar (jMaciejewskill2003l . l2004j ; lAnn fc Thakurll2005l ). Thus, it is quite plausible that the 
central mass concentration 7 can affect the gas inflow into central hundred parsec of models 
M1-M4 with c s = 10 km s _1 , since Figure 3 shows that the central mass concentrations of 
7 = and 0.5 make physical conditions in which it is plausible for the symmetric two-armed 
nuclear spirals to take the form of hydrodynamic shocks. This inspires us to examine the 
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effect of the central mass concentration 7 on the gas inflow into the central hundred parsec 
of models M1-M4 in detail. In order to do this, the time evolution of the fraction of gas 
particles accumulating inside the radial zone of 0.3 kpc is shown in Figure 4 for models 
M1-M4. The gas accumulated inside the radial zone of 0.3 represents the gas inflow deep 
inside the ILR, since the OILR in model Ml and the ILRs in models M2-M4 are located 
around at 1.3 kpc and 1.5 kpc, respectively. As can be seen in Figure 4, models M1-M4 show 
appreciable gas inflow into the central hundred parsec throughout the evolution. Although 
there are little differences in the gas inflow among the models in the early phase of evolution, 
model M4 with 7 = 1.5 shows slightly larger gas inflow up to ~ 6rt, ar . However, at the later 
evolution times (i.e., after ~ 6r;, ar ), the remarkable differences in the gas inflow are found 
in the sense that the gas inflow increases as the central mass concentration decreases from 
7 = 1.5 to 0. As a result, the highest gas inflow rate is found for model Ml with 7 = 
whose mass is less centrally concentrated than models M2-M4 that assume 7 = 0.5, 1 and 
1.5, respectively. This correlation between the amount of gas inflow close to the center and 
the central mass concentration still exists when the models with each of the remaining values 
of the gas sound speed (i.e., c, — 15 and 20 km s" 1 ) a re considered. However, similar to 



Patsis fc Athanassoulal (120001 ) and lAnn fc Thakurl (120051 ) , the gas inflow rate increases when 



the gas sound speed increases from c s = 10 to 20 km s . 



4. DISCUSSION 

4.1. Dependence of Gas Flow Morphology on Central Mass Concentration 

and Gas Sound Speed 



Using the SPH method, lEnglmaier fc Gerhard! (119971 ) have found that when the gas 



in barred galaxies is modelled as an isothermal fluid, the gas flow morphology depends on 
the shape of the underlying potential and the sound speed. They noticed that the off-axis 
principal shocks with a nuclear ring developed by the gas flow at low sound speed turn into 
the on-axis ones with a straight inflow to the center when the sound speed is high. However, 
the resolutions of their simulations were not high enough to analyze the nuclear morphology 
in detail. Our present SPH simulations confirm their findings that not only the underlying 
potential, which can be greatly affected by the central mass concentration, but also the gas 
sound speed influence the morphology of the gas flow. 

Because of the better resolution of our SPH simulations, we analyze the gas flows in the 
central kiloparsecs of barred galaxies, where the morphology of nuclear spirals is found to be 
dependent on the central mass concentration 7 and the gas sound speed c s (see § 3.1). The 
openness of the nuclear spiral arms increases as the central mass concentration decreases. 
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Increasing the gas sound speed has a comparable effect on the openness of nuclear spirals to 
that of decreasing the central mass concentration. This implies that the nuclear spiral arms 
also open out when the gas sound speed increases. As a result, the less centrally concentrated 
model with high gas sound speed (i.e., model M8 with 7 = and c s = 20 km s _1 ) shows more 
loosely wound nuclear spiral whose openness appears to be larger than those of the other 
selected models. This can be justified by the fact that the innermost parts of the nuclear 
spiral formed in model M8 extend all the way to the galactic center with winding angle of 
~ 3tt (see § 3.1 and Fig. 2). A similar effect of the central mass concentration and the 



gas so und speed on the morphology of nuclear spirals was found by lEnglmaier Shlosman 
( I2OOOI ). although th ey used the two-dimensio nal grid-based hydrodynamics code ZEUS-2D. 
In addition to this, iMaciejewski et all (120021 ). who use the grid-based hydrodynamics code, 
have also found that the nuclear spiral opens out in the sense that its innermost parts reach 
very close to the galactic center when the gas sound increases from c s = 5 to 20 km s^ 1 . 



Ann fc Thakurl (120051 ). who use the SPH code, have found that the models with lower 
gas sound speeds allow a larger diameter of the nuclear disk than those with high gas sound 
speeds. Here we have also found a similar effect of the gas sound speed on the diameter of 
the nuclear disk (see § 3.1 and Fig. 2). Beside this, we have noticed that the diameter of 
the nuclear disk increases as the central mass concentration 7 increases. A similar effect of 
the central mass c oncentration and the gas sound speed on the diameter of the nuclear disk 
has been found by lEnglmaier fc Shlosman! (120001 ). Thus, it is apparent that increasing the 
central mass concentration has a similar effect on the diameter of the nuclear disk to that of 
decreasing the gas sound speed. 



Furthermore, it was also found by lAnn fc Thakurl (120051 ) that the nuclear regions of the 
modeled gaseous disks evolve to trailing nuclear spirals whose morphologies depend on the 
gas sound speed in the presence of a central SMBH that is massive enough to remove the 
IILR, whereas the lower gas sound speed model (c s = 10 km s -1 ) that has the existence 
of the IILR due to the absence of a central SMBH evolves to a leading nuclear spiral. In 
general, the leading nuclear spir als are generated at the IILR because of the positive torque 
of the bar (jCombes et al.ll2002l ). while the trailing nuclear spirals at the OILR (or at the 
ILR when there is no IILR) are made by orbital switchi ng from the xj -orbi t s that align 



with the bar axis to the x 2 -orbits perpendicular t o the bar (lAthanassoulalll992l ; l\Vadalll994 



Knapen et al.lll995l ; lEnglmaier fc Shlosmanll2000l ). We have not seen here any such existence 
of the leading nuclear spirals, since the masses of models with 7 = are sufficiently centrally 
concentrated to weaken the IILR by placing its location at the negligible distance from the 
galactic center (see Fig. 1). This makes physical conditions plausible for the development 
of trailing nuclear spirals in the models with 7 = whose innermost parts can reach the 
galactic center when the gas sound speed is high (see Fig. 2). A similar dependence of the 
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morphology of trailing nuclear spirals on the gas sound speed was found for the models with 
7 > 0.5 that have only a single ILR due to the increase of central mass concentration from 
7 = 0. 

The relative bar strength is important in galaxy morphological studied because phenom- 
ena such as the gas inflow can be affected in various ways to the eff ectiveness with which 
the bar potential influences the motion o f gas in a galactic disk fee. . ISe 



loow &: Wilkinson 



19931 : iButa & Combeslll996l : lKnapenlll999f ). As mentioned in Block et all (1200 If ), the relative 
bar strength defined by the maximum of the tangential-to-radial force ratio (F Tan / F Rad ) max 
increases when the bulge becomes weak. We have found here that the gas inflow into the 
central hundred parsec increases as the central mass concentration decreases from 7 = 1.5 to 
(see § 3.3 and Fig. 4). This correlation between the amount of gas inflow close to the center 
and the central mass concentration can be understood if we consider the dependence of the 
relative bar strength (Fran/ ' F^ad] 'max on the strength of the bulge. Since the relative bar 
strength (Fran/ Fn a d) max increases as the bulge becomes weak with the decrease in central 
mass concentration 7, the gravitational torques of the bars would be stronger in the models 
with 7 = than in the models with 7 = 1.5. This provides favorable conditions for the 
gas particles in the models with 7 = to lose their angular momentum much faster than 
those in the models with 7 = 1.5. Thus, the highest gas inflow rate is found for model Ml 
with 7 = (see § 3.3 and Fig. 4), since the gas particles move inward because of the loss of 
angular momentum. Apart from the relative bar strength, the gas inflow close to the center 
can also be understood by the ILR that exists only if the central mass concentration is large 
enou gh (lEnglmaier &: Gerhardl 119971 ; iPatsis fc Athanassoulal I200CH : lEnglmaier &: Shlosman 
2000). The ILR is defined as the outer limit to which the X2-orbits can extend (IMaciejewski 
20041 ). If th e ILR is absent, the X2-orbits do not ex ist, and strong gas inflow to the center can 
take place (lAthanassoulal Il992l ; iPiner et al.lll995l ). As can be seen in Figure 2, decreasing 
the central mass concentration from 7 = 1.5 to corresponds to the decrease in the strength 
of the ILR in the sense that the diame ter of the nuclear disk, whic h represents the region 
occupied by the o^-orbits, decreases (cf. lEnglmaier fc Shlosmanll2000l ). This causes the high- 
est gas inflo w rate in the less central l y con centrated model (eg., model Ml with 7 = 0), as 
suggested by lEnglmaier fc Shlosmanl (120001 ) . 



Using the SPH method, IPatsis fc Athanassoulal (120001 ) and lAnn fc Thakurl (120051 ) have 



found that the gas inflow inside the radial zone near the central regions increases with the gas 



(Enelmaier & Gerhard 


1997 


) 


2000; 


Maciejewski et al 


2002 



9970 or the grid-based hy drodynamics code (lEnglmaier fc Shlosman 



20021 ; iMaciejewskil 120041 ) , also indicate that the gas inflow close to 



the center increases with the gas sound speed. Although we have only presented here the 
effect of the central mass concentration on the gas inflow, a closer look at Figure 2 suggests 
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a similar correlation between the amount of gas inflow close to the center and the gas sound 
speed to t hat found in previous simula tions. This can be understo o d if w e consider the 
findings of iPatsis &; Athanassoulal (120001 ) and lEnglmaier fc Shlosmanl (120001 ) that while the 
gas particles in the low sound speed mediums easily follow the x 2 -orbits, the gas particles 
in the high sound speed mediums face more difficulty doing so. Because of this reason, 
Figure 2 shows that the diameter of the nuclear disk, which represents the region occupied 
by the a^-orbits, decreases when the gas sound speed increases from c s = 10 to 20 km s _1 
(cf. lEnglmaier fc Shlosmanl l2000l ) . resulting in a large amount of gas inflow deep inside the 
ILR for the models with high gas sound speed (c s = 20 km s~ x ). Thus, the above described 
correlations make physical conditions in which it is plausible for the innermost parts of 
nuclear spiral to reach very close to the center when the less centrally concentrated model 
with high gas sound speed is considered (see Fig. 2). 



4.2. Mechanism Responsible for the Formation of Nuclear Spirals 



hydrodynamic shocks (IMaciejewski et al. 



The gas density w aves (lEnglmaier fe Shlosmanl l200d; IMaciejewski et al.ll2002l ) and the 



2002 



Ann fc Thakurl 120051 ) are two promising 



mechanisms proposed for t he formation of the symmetric two- armed nuclear spira l s in b arred 
galaxies. As mentioned in lEnglmaier fc Shlosman (hoooh and iMaciejewski et all ( 120021 ) , the 
linear theory of gas density wave postulates the following: (1) the spiral pattern opens out as 
the gas sound speed increases; (2) the spiral pattern does not extend beyond the IILR when 
two ILRs are present; (3) when one ILR is present, the spiral pattern can exist between 
the center and the ILR. Since the gas density waves cannot propagate in regions where 
shocks are strong enough to drive a nonlinear response of gas flows, it can be applied to 
low-ampli tude spirals where no shocks exist and the gas respon se can be considered a lin- 
ear wave (jEnglmaier Sz Shlosman 2000 ; Maciejewski et al. 2002 ). However, it is worthwhile 
to mention that the gas densit y wave itself is triggered by t he spiral shocks driven by the 
gravitational torque of the bar (lEnglmaier fc Shlosmanl l2000l ). 



The grid-based hydrodynamic simulations of IMaciejewski et al.1 (120021 ) show that the 
gas sound speed plays a significant role to decide whether the nuclear spirals are formed by 
gas density waves or by hydrodynamic shocks. In their low gas sound speed run (c s = 5 km 
s -1 ), the nuclear spiral is weak and rather tightly wound that creates at most a factor of 
2 density enhancement over the interarm region. Moreover, it terminates just outside the 
IILR, and its shock strength, as measured by div 2 v for div v < 0, is 10 times smaller than 
that of the principal shock. This suggests that the nuclear spiral in their low gas sound speed 
run (c s = 5 km s _1 ) can be well understood in terms of the gas density waves. On the other 
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hand, the strong shock nature of nuclear spiral is found by lMaciejewski et all ( 120021 ) in their 
high gas sound speed run (c s = 20 km s _1 ). They have found that the strong loosely wound 
nuclear spiral is a direct continuation of the principal shock along the leading edge of the 
bar, which generally agrees with the linear t heory of gas density wave tha t postulates a open 
nuclear spiral for a high gas sound speed (cf. lEnglmaier fc Shlosmanll2000l ). Nevertheless, the 
nuclear spiral reaches close to the nucleus well beyond the IILR where the gas density waves 
cannot penetrate. In addition to this, the shock strength (i.e., div 2 v for div v < 0) and the 
arm/interarm density contrast of the nuclear spiral are comparable to those of the principal 
shocks. This indicates that the nuclear spiral in their high gas sound speed run (c s = 20 km 
s _1 ) is definitely well beyond the linear regime explored in the density wave theory where 
sound waves cannot propagate in regions for w hich shocks are stron g. Our previous SPH 
simulations in the high sound speed mediums (lAnn &: Thakurl 120051 ) also show the strong 
shock nature of nuclear spirals. In these simulations, we have found that the symmetric 
two-armed nuclear spirals are formed by hydrodynamic shocks in the presence of a central 
SMBH that can remove IILR when the gas sound speed is high enough to drive a large 
amount of gas inflow deep inside the ILR. 



Like the role of gas sound speed as described above (cf. lMaciejewski et al.ll2002l ; lAnn fc Thakur 



20051 ). we have found here that the central mass concentration also play a significant role in 



the mechanism responsible for the formation of the symmetric two-armed nuclear spirals in 
barred galaxies. Our present SPH simulations show that the symmetric two-armed nuclear 
spiral arms open out as the gas sound speed increases (see § 3.1 and § 4.1), which gener- 
ally agrees with the predictions of the linear theory of gas density wave. Nevertheless, the 
symmetric two-armed nuclear spirals of the models with 7 = and 0.5 cannot be explained 
by gas density waves, since they are supported by the gas particles that are shocked by 
supersonic flow (see § 3.2 and Fig. 3). On the other hand, since the symmetric two-armed 
nuclear spirals of the models with 7 = 1 and 1.5 cannot be explained by the shocked gas 
particles (see § 3.2 and Fig. 3) and their arms open out when the gas sound speed in- 
creases, they can be well understood in terms of the gas density waves. This suggests that 
the symmetric two-armed nuclear spirals formed by hydrodynamic shocks turn into those 
explained by gas density waves when the central mass concentration increases from 7 = 
to 1.5. Since the relative strength of bars in the models with 7 = is stronger than that in 
the models with 7 = 1.5 (see § 4.1 ), our results appea r to confirm the findings of the grid- 
based hydrodynamic simulations of iMaciejewskil (12004 ) that the nuclear spirals developed in 
the models with strong bars have the nature of shocks in gas, whereas those formed in the 
mod els with weak ovals can be explained in terms of the gas density waves. In addition to 
this, IMaciejewskil (120041 ) seems to support our findings that the nuclear spirals formed by 
hydrodynamic shocks trigger a large amount of gas inflow than those explained in terms of 
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the gas density waves. 



Although our present SPH simulations and the grid-based hydrodynamic simulations of 
Englmaier fe Shlosmanl (120001 ) show a similar effect of the central mass concentration and 
the gas sound speed on the mor phology of nuclear spir a ls (see § 4.1), the shock nature of 
nuclear spiral was not found by lEnglmaier &: Shlosmanl (120001 ) . They have found the low- 
amplitude nuclear spirals, which can be explained by gas density waves because the spiral 
arms open out as the gas sound speed increases. On the other hand, we have found here 
that the nuclear spirals can be formed either by hydrodynamic shocks or by gas density 
waves, depending on the value of central mass concentration 7. It appears that the different 
assumptions of the mass models along with the adoption of the different hydrodynamics 
codes m ight play some role in cr e ating this discrepancy between our present simulations and 
those of lEnglmaier fe Shlosmanl (|2000| ). 



4.3. Effect of Number of Particles 

After performing the SPH simulations with 2 x 10 4 and 1 x 10 5 SPH particles for 
the response of the g a seous disk to the imposition of nonaxisymmetric bar potentials, 



Patsis fc Athanassoulal (120001 ) claimed that the morphology of gas inflow does not depend 
on the number of SPH particles for N g larger than 2 x 10 4 . This encourages us to check 
whether the number of SPH particles causes any major changes in our results discussed in the 
previous sections. Thus, we have repeated our simulations for models M1-M10 with 1 x 10 6 
SPH particles that give the initial resolution of ~ 20 pc. Since we are mainly interested 
to examine whether the number of SPH particles affects our results discussed in §§ 3.2 and 
4.2 regarding the role of the central mass concentration 7 in the mechanism responsible for 
the formation of the nuclear spirals in barred galaxies, we have presented here the nuclear 
distributions of supersonically shocked gas particles (i.e., — /iV ■ v > c s ) for models M1-M4 
with 1 x 10 6 SPH particles at the evolution time of 10Tb ar hi Figure 5. Similar to Figure 
3 that depicts the nuclear distributions of supersonically shocked gas particles for models 
M1-M4 with 2 x 10 4 SPH particles at the evolution time of 10r 6ar , Fi gure 5 also shows that 
the nuclear spirals formed in models Ml with 7 = and M2 with 7 = 0.5 have the nature 
of hydrodynamic shocks in gas, whereas those developed in models M3 with 7 = 1 and M4 
with 7 = 1.5 are not supported by the shocked gas particles. Although it is not presented 
in this paper, the results of our simulations, which employ 1 x 10 6 SPH particles, are sim- 
ilar to those obtained in the previous sections using 2 x 10 4 SPH particles for the effect of 
the central mass concentration 7 and the gas sound speed c s on the morphology of nuclear 
spirals (see § § 3.1 and 4.1), as well as on the gas inflow close to the center (see § § 3.3 and 
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4.1). This implies that similar to our simulations with 2 x 10 4 SPH particles (see § § 3.1 
and 4.1), the nuclear spiral arms formed in all the models with 1 x 10 6 SPH particles also 
open out as the gas sound speed c s increases. Because of this fact, our simulations with 
1 x 10 6 SPH particles suggest that the nuclear spirals of models M3 with 7 = 1 and M4 
with 7 = 1.5 can be explained in terms of the gas density waves (cf. lEnglmaier fc Shlosman 



20001 ) . which is consistent with the formation mechanism proposed in § 4.2 for the nuclear 
spirals of these models with 2 x 10 4 SPH particles. As a result, we do not see any changes 
in our findings regarding the effect of the central mass concentration 7 on the mechanism 
responsible for the formation of the nuclear spirals in barred galaxies when the number of 
SPH particles increases from 2 x 10 4 to 1 x 10 6 . From this, we can conclude that similar to 



Patsis fc Athanassoulal (120001 ). our results do not depend on the number of SPH particles for 



N g larger than 2 x 10 4 , and thus that this number of SPH particles is sufficient to examine 
the shock nature of nuclear spirals in barred galaxies. 



5. CONCLUSIONS 

We have used the SPH method to find the response of the nuclear region of a gaseous 
disk to the imposition of nonaxisymmetric bar potentials in order to understand the effect 
of the central mass concentration 7 on the mechanism responsible for the formation of the 
symmetric two-armed nuclear spirals in barred galaxies. We have found that the central 
mass concentration 7 plays a significant role in shaping the gravitational potential in the 
central kiloparsec to weaken or to remove the IILR that prevents the gas inflow close to the 
nucleus. It was also noticed that the gas flow driven by the gravitational torque of the bar 
leads to the formation of symmetric two-armed nuclear spiral whose morphology depends on 
the central mass concentration 7 and the gas sound speed c s in the sense that the nuclear 
spiral arms open out when either the central mass concentration decreases or the gas sound 
speed increases. Thus, the less centrally concentrated model with high gas sound speed 
shows more loosely wound nuclear spiral whose innermost parts reach the galactic center in 
a curling manner. A largest nuclear disk is found for the more centrally concentrated model 
with low gas sound speed. 

Furthermore, the central mass concentration 7 is found to play a crucial role to decide 
whether the symmetric two-armed nuclear spirals are formed by hydrodynamic shocks or by 
gas density waves. The symmetric two-armed nuclear spirals are formed by hydrodynamic 
spiral shocks caused by the gravitational torque of the bar for the models with 7 = and 
0.5, whereas they are explained in terms of the gas density waves when the models with 
7 = 1 and 1.5 are considered. This suggests that as the central mass concentration increases 



-22 - 



from 7 = to 1.5, the mechanism responsible for the formation of the symmetric two-armed 
nuclear spirals in barred galaxies shifts from the hydrodynamic shocks to the gas density 
waves. The symmetric two-armed nuclear spirals formed by hydrodynamic shocks in the 
models with 7 = drive a large amount of gas inflow than those explained by gas density 
waves in the models with 7 = 1.5. Finally, we conclude that 2 x 10 4 SPH particles give the 
same results as 1 x 10 6 SPH particles, and thus are adequate to examine the shock nature 
of nuclear spirals in barred galaxies. 
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Fig. 1. — Rotation curves and angular frequencies of model galaxies. Because of the non- 
axisymmetric nature of the bar, its contribution is included here after averaging it axisym- 
metrically. Left: Model Ml with 7 = 0; right: model M3 with 7 = 1. The horizontal 
dot-dashed lines in the bottom panels represent the bar pattern speed of Q p = 42.8 km s _1 
kpc -1 . 



-26- 



7=0.0 



-1 
1 

7=0.5 

-1 
1 

7=1.0 

-1 
1 

7=1.5 



t~ 1 — 1 — r 

"10 



10 



10 



Ml" 



H 1 1 h 



M2~ 



H 1 H — I H 



M3 



H h — ■)— I- 



10 



M4 



_i i j 1 



15 



M5 



H h— H h— 4- 



15 



M6~ 



I I — I h 



15 



M7 



_j — 1 1 1_ 



III 



20 



20 



20 



M8~ 



I I 1 h 



M9~ 



-+— I h 



M10 



- -1 



1 -1 



Fig. 2. — Nuclear features of the gaseous disks in models M1-M10 at the evolution time of 
10T{, ar . The bar rotation period Tf, ar is ~ 1.44 x 10 s yr. The models are indicated in the top 
right corner and the gas sound speed in units of kilometers per second is given in the top 
left corner of each panel. On the left side of each row, its corresponding value of the central 
mass concentration 7 is given. The bar lies horizontally, and the box size is 2.4 kpc in one 
dimension. 
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Fig. 3. — Nuclear distribution of the gas particles that are shocked by supersonic flow (i.e., 
—hV ■ v > c s ) for models M1-M4 at the evolution time of 10rb ar . Red: model Ml with 
7 = 0; blue: model M2 with 7 = 0.5; green: model M3 with 7 = 1; magenta: model M4 
with 7 = 1.5. The gas sound speed is c s = 10 km s -1 . The number in the top left corner of 
the panel is the evolution time in units of Tb ar . The bar lies horizontally, and the box size is 
2.4 kpc in one dimension. 
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Fig. 4. — The inflow of the gas particle into the central hundred parsec of models M1-M4. 
The time evolution of the fraction of gas particles accumulated inside a radius of 0.3 kpc is 
shown. The time is given in units of r& ar . The gas sound speed is c s = 10 km s . The line 
colors corresponding to each model are given in the top left corner of the panel. 
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Fig. 5. — Same as Fig. 3, but for 1 x 10 6 SPH particles. 



